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Abstract 



We consider the Wigner ensemble of Hermitian n-dimensional random ma- 
trices W}j^ = Wij /\/n and study the asymptotic behavior of the expression 

K n (s', s") = E jlr (W in) ) 2s ' Tr (W (n) ) 2s " \ - ETr (W (n) ) 2s ' ETr (W {n) ) 2a " in 

the limit s',s",n — > oo such that s' — s' n and s" = s'^ are the values of the 
order n 2//3 . Assuming that the random variables {wij} are of the symmetric 
probability distribution such that its all moments V^k are of the sub-gaussian 
form, we prove that the limit of K n (s' n , s„) exists and does not depend on the 
particular values of V^h, k > 2. 

The proof is based on the combination of the arguments by Ya. Sinai and 
A. Soshnikov with the detailed study of a moment analog of the Inverse Partic- 
ipation Ratio of the Gaussian Unitary Invariant Ensemble of random matrices 
(GUE). 

Running title: Correlation Function of High Moments 

1 Introduction 

The famous semi-circle law is proved by E. Wigner in the fifties by the moment 
method [13]. Since that, the moment method is widely used in the spectral theory 
of large random matrices. In particular, it has been employed in the studies of the 
asymptotic behavior of the spectral norm [1, 4, 5], where one considers the moments 
of the order that grows at the same time as the matrix dimension tends to infinity. 

Later, the moment method has been developed to describe the properties of the 
eigenvalue distribution on the local scale [11] that gives much more detailed infor- 
mation than that of the semi-circle law known as the global (or integral) one. In 
the subsequent paper [12], the so-called microscopic scale at the edge of the limiting 
spectra is reached in the frameworks of the proof of the universality conjecture of 
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the properties of the local eigenvalue distribution. In [12], a spacious program of 
the studies of the high moments of the Wigner ensemble as well as their correlation 
functions (or cumulants) has been proposed and strong method has been developed. 

In our previous two papers, first [8] and then [7], we followed this way and con- 
sidered the asymptotic behavior of high moments of large Wigner random matrices. 
In the present paper we continue the work in this area and pass to the studies of the 
second correlation function (or the second cumulant) of these moments. 

2 Main results 

Let {Xij,Yij, 1 < i < j} be family of jointly independent real random variables 
determined on the same probability space and let E denote the mathematical expec- 
tation with respect the corresponding probability measure. Assuming that 

EXij = EYij = 0, (2.1) 

and 

EX£ = (l + J tf )/8 and EY? = (I - 5^/8 (2.2) 

for all i and j, where Sij denotes the Kronecker 5- function, we determine complex 
random n x n matrices with the elements 

^-Mx»-%**iii (2 ' 3) 

and say that the family {W( n )} represents the Wigner ensemble of Hcrmitian random 
matrices [13]. Below we will omit the superscripts n when no confusion can arise. 
We are interested in the asymptotic behavior of the covariance of the traces 

n 
i=l 

in the limit when s and n tend to infinity. Our main result is as follows. 

Theorem 2.1. Let random variables Xij and Yij have the same symmetric 
probability distribution such that, in addition to (2.2), the moments E(Xij) 2k = V2k 
exist for all k > 2 and Vik < (cfc) fe for some c > 0. Then there exists xo > such 
that the limit 

K(x', X") = lim (E |Tr W 2s ' Tr W 2s " j - ETr W 2s ' ETr W 2s "^j , (2.4) 

where 

S ' = 4=L(x'n 2 ) 1/3 J, *" = <=L(x"n 2 ) 1/3 J, (2.5) 

for any given positive \' an d x" ^ ess than xo, exists and does not depend on the 
particular values ofV2k, k>2. 

The proof of Theorem 2.1 is performed consists of two parts. On the first step 
we reduce the study of the covariance 

K n {s', s") = E jlr W 2s ' Tr TT 2s " j - ETr W 2s ' ETr W 2s " (2.6) 
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to the study of a part of an expression that can be regarded as a moment analog 
of the Inverse Participation Ratio of the Wigner ensemble (2.3). We do this with 
the help of the arguments used in paper [11] combined with the results on the upper 
bound of the moments of the Wigner ensemble obtained in [8] . Let us note that the 
condition of the sub-gaussian moments Vik < (ck) k imposed in Theorem 2.1 can be 
essentially relaxed [7]. We keep it here to make the main stress on the second part 
of the proof of Theorem 2.1. 

The second stage contains a detailed study of the moment analog of a Green 
function version of the Inverse Participation Ratio considered in the case of the 
well-known Gaussian Unitary Invariant Ensemble (GUE). This ensemble represents 
a restriction of the Wigner ensemble (2.3) with (2.1) and (2.2) to the case of the joint 
Gaussian (normal) distribution of the random matrix elements [9] . On this stage we 
use the method of non-asymptotic estimates of the moments of GUE developed in 

Everywhere below, we will say for the sake of simplicity that K n (s',s") (2.6) 
represents the correlation function of the moments of random matrices W^ n \ 

3 Correlation function and moment analog of IPR 

The study of the correlation functions of high moments of Wigner random matrices 
has been started in [11]. The reasoning of this paper is based on the natural repre- 
sentation of K n (s', s") (2.6) as a weighted sum over the set of pairs of closed paths 
of 2s' and 2s" steps 

where /« = (#>, . . . , i™^), 1^ - (i< 2) , . . . , i^-i,^). and *f e {1, 2, . . . ,n} 
for k = 1 , 2 and all j ; namely 

K n {s',s")= ]T U n (I^), (3.1) 

' (1 ' 2) ^;?L» 

where Il n is a weight of the pair j( 1,2 ) given by the mathematical expectation of the 
product of random variables Wij/y/n determined by I^ 1,2 \ 

Clearly, a non-zero contribution to the right-hand side of (3.1) is given by a subset 
of path pairs such that 1^ and 1^ have at least one step in common. This means 
that there exist two instants of time t' and t" such that 

= #> and = #> +1 } or = ^ +1 and ^ = #>} ( 3 .2) 

(here and below we identify i^J, with and i^}„ with i^). In these two cases, 
we say that the step (t',t' + 1) of 1^ is passed by 1^ in the direct and the inverse 
orientations, respectively. In what follows, we consider t' and t" to be the first instants 
of time, i.e. the minimal ones, with the property (3.2) to hold. The number of times 
that the step (t', t' + 1) is seen in the path pair is referred as to the multiplicity \i of 
the step. 

In paper [11], the main attention is paid to the subset Z^ls"^) °f P a th pairs 
such that 1^ and 1^ have at least one step of multiplicity \i = 2 in common. This 
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kind of path pairs is referred to as to the simply correlated pairs [11]. Since 

are Hermitian, this step in common is passed by 1^ in the inverse orientation with 

respect to 

By using the arguments of [11], one can prove that the sum over the set of simply 
correlated paris 

t n (s',s")= ]T IIn(/ (1,2) ) (3-3) 

/(1 ' 2)ei ^ 2 L» 

remains bounded in the limit (2.5) asn^oo and that the contribution of the path 
pairs such that there exists at least one step of multiplicity [i > 4 vanishes, no matter 
how many times this step is seen in 

Indeed, regarding a simply correlated pair one can remove the steps 

(t',t' + 1) and (t",t" + 1) and consider the collection of remaining steps as a closed 
path of 2s' + 2s" - 2 steps 

T _ (aW ,-(!) ,-( 2 ) ;( 2 ) J 2 ) ,-(2) •(!) ,-Uh 

- i 2s'+2s"-2 — \ l o ' ' • • ' V ' f t"+2' ■••>'() ' f l ' ' • - ' f t" ' V+2' ■ • • ) 4 / 

such that 

n(^;i») = ^n(/ 2s , +2s „_ 2 ). (3.4) 

It is natural to say that the two-steps reduction procedure (1^,1^) — > J 2s '+2 S "-2 
is performed here. 

To get a non-zero weight n(/ 2s / +2s //_ 2 ), one has to consider /2s'+2s"-2 to be an 
even closed path [11]. Since 

J2 n(/ 2s , +2s „_ 2 ) = E {Tr ( W (")) 2s '+ 2s "- 2 } , 

^2s' +2s" -2 

one can use the result about the universal upper bound for the averaged moments 
nM 2 ( " ) +2s „_ 2 = E jTr(VF(")) 2s '+2s"-2j [12] and prove the existence of their limit. 

Let us stress that the proof of this estimate given in [12] should be completed and 
modified (see [7] or [8]). To show that T, n (s',s") remains bounded, it suffices to 
estimate the number of simply correlated path pairs i^ 1 ' 2 ) that can be reconstructed 
from an even closed path l2s'+2s" -2- Omitting the details, one can say that this 
number is proved in [11] to be proportional to (2s' + 2s" — 2) 3 / 2 that is compensated 
by the factor l/(4n) of (3.4) and this completes the argument. 

As for the path pairs that are not simply correlated ones, it is claimed in [11] 
and subsequent papers that the study of the corresponding sums can be also trivially 

(n) 

reduced to the study of Af 2s / +2 s"-2- Indeed, using the arguments described above, 
it is not to hard to show that the contribution to the right-hand side of (3.1) that 
comes from the path pairs such that and /( 2 ) have at least one step of multiplicity 
fi > 6 in common, vanishes in the limit (2.5) as n — > oo. The same concerns the path 
pairs that have one step of multiplicity fj, = 4 in common and have some other step 
in common passed 4 times or more. 

However, it remains one more case of non-simply correlated pairs that seems to 
be not so easy to treat. Consider a path pair such that and 1^ have a step 
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of multiplicity 4 in common and do not have any other common step. Then the 
two-steps reduction procedure described above leads to a path /2s'+2s"-2 that, in 
particular, can be free from the steps of multiplicity fj, > 4. Then it is not clear how 
to prove the fact of vanishing contribution of the corresponding sum E„(s',s") by 
the use of the arguments of [11] only. Up to our knowledge, no rigorous study of this 
question has been published or reported. In the present paper we give a proof based 
on the method of [6]. 

If the common step is passed by 1^ twice in opposite directions, then we get a 
sub-sum of S„(s', s") given by expression 

^V,*") = § Y.Y.vi^^^} V{\W^-^W^} . (3.5) 

Here the symbol [ • ] * indicates the fact that the corresponding product does not 
contain the steps of multiplicity 4 or greater. Therefore we can bound the right-hand 
side of the previous inequality by the expression 

Ms',s") = ^ E E E — { A %r A li,} Eoue {A%- v Al i2 } , (3.6) 

where A = are the random matrices of the Gaussian Unitary Invariant Ensemble 
(GUE) and E GUE denotes the corresponding mathematical expectation [9]. Certainly, 
we assume that the matrix elements A^ have the mean zero and the variance (4n) _1 
(cf. (2.2)). We used here the fact that all terms of the right-hand side of (3.6) are 
positive. 

If the common step is passed by twice in the same direction, we get a sub-sum 

U,V i\,l2 

bounded by 

wo = § EE E OUE {4 s ;r^4 e gue {4*;;-^ 2 }, (3.7) 

In the proof, we will see that S n (s' , s") — o(R n (s' , s")) in the limit (2.5) as n — > oo. 
If the common step is passed by one or three times, the corresponding sub-sums 
are zero. 

Let us point out that the expression (3.6) resembles the expression 

n 

J2^{G xx (zi) G xx (z 2 )} (3.8) 

x=l 

with G{z) = (A — zl) , z\ = A + 'vq, z 2 = A — vq, r\ > known in theoretical 
physics as a version of the Inverse Participation Ratio. This quantity that reflects the 
(dc)localization properties of the eigenvectors of random operators naturally appears 
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also in the studies of the spectral properties of random matrices [2, 10]. Recently, 
it has been studied in the frameworks of the proof of the universality of the bulk 
spectral distribution of Wigner random matrices [3] on the microscopic scale 1/n. 

It is natural to consider (3.6) as a natural analog of (3.8). This become even 
clearer when one considers the diagonal part of (3.6) (see Section 5). Let us note 
that in the studies of the edge of the spectrum of Wigner ensemble, the asymptotic 
regime (2.5) corresponds to the microscopic scale n~ 2 / 3 . 

4 Estimates of moments of GUE 

Let us consider the Gaussian Unitary Invariant Ensemble of random matrices (GUE) 
that is given by a family of random Hermitian matrices with elements 

U (n) ) =^a xyi x,y = l,...,n (4.1) 



such that the law of A^ n > has a density proportional to exp{— 2nTr (A^ n ') 2 } [9]. Then 
(2.1) and (2.2) hold. For p non-negative integer, we denote 

n 1 

L™ = -^(AWy xx = -Tr(AMy>, (4.2) 

x=l 

and consider the averaged moments Mp"-* = ELp™' , where E denotes the mathematical 
expectation with respect to the probability measure generated by {A^™'}. We also 
denote Up n \x) = ~Ei(A^)p. x . Everywhere below, we will omit the superscripts n 
when no confusion can arise. 

We represent R n {s' ', s") (3.6) as a sum of four terms 

Rn(s',s") = ^ £J#V>A (4-3) 
fc=l 

such that 

R£\s',s")= £ ( Yl U ai (x)U Pl (y)) ( E U a2 (x)U 02 (y) 

x,y=l \a 1 +P 1 =2s> J \a 2 +/3 2 =2s" 

W (*',*")= E ( E U ai (x)U Pl (y)\ I Yl E{(A-r(^)°}], 

x,y=l \ai+/3i=2s' J \a 2 +fe=2s" 

4 3) (*V) = E ( E E{(A-)°(0°}] ( E ^O^fo) 

*,3/=l \ai+/3i=2s' / \q 2 +/3 2 =2s" 

and 

rW(s',s")=y[ E E{(A-r«r})( ^ E{(A-r(^r}|, 

x,y=l \ai+ j 8i=2s' / \a 2 +/3 2 =2s" 
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where we denoted by X° the centered random variable X° = X — EX. 



In the present section we study the values U2 S (x). We do this mostly in the 
frameworks of the method of recurrent non-asymptotic estimates developed in [6] 
with respect to the moments M^" of GUE. The statement we prove generalizes the 
results of [6] and can be formulated as follows. 

Theorem 4.1. Given any constant h > 1/16, there exists < k < 12 — 3/(4/i) 
such that the estimate 

( s(s 2 -l)\ 
sup U 2s {x) < (l + h V > \m 8 (4.4) 

x=l,...,n V " I 

holds for all values of integer positive s and n satisfying condition s 3 /n 2 < k, where 
m s are the moments of the corresponding to (4-1) semi-circle distribution [13] 

1 {2s\ 1 (2sV 

m ^¥^TT){ s )-^-^diyr -0,1,2,... (4.5) 

Regarding the generating function tp{r) = X)s*lo msTS > one can rewr ^ e (4-4) in the 
form 

[ hr 2 



sup U 2s (x) < 

x — 1,. ...n 



:(l-r)«/» ■ (4 ' 6) 
where [/(t)] s = f s denotes the coefficients of the corresponding generating function 

The proof of Theorem 4.1 is based on the analysis of recurrent relations for U2 S (x) 
and related variables 

D$(x)= E \V{(A a r xx K 2 ---L° ar }l «i>l- 

ai + ...+a r -2s 

In fact, we need to consider more general than U2 S (x) variable 

n 

U 2s (x, y) = ~E(A 2s ) xy = E E {A xt (A 2s ~ 1 ) ty } . 
t=i 

Applying to the last mathematical expectation the integration by parts formula (see 
Section 7 for the details), we get equality 

n 1 to 2s — 2 

t=l t=l j=0 



= ^ M ^^(^)+ E E{(A^)°,i° 2 }- (4.7) 

j— ai+ft2- 2s — 2 

Introducing variable 

Dt\x,y)= E |E{(^)°,i° 2 ---i° r }|, 

ai+...+a r =2s 
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we deduce from (4.7) the main inequality 



U 2s (x,y)<-(M*U(x,y)) 2s _ 2 +D 2 2 J_ 2 (x,y), s>l (4.8) 



with the initial condition Uq{x, y) — 5 xy , = 0. 
Regarding the expectation 

n 

E { (A a > ) xv [Ll 2 ■ ■ ■ L° a r ] ° } = E {A x t {A ai ) ty [L° a 2 • • • L° r ] ° } 

t=i 

and using again the integration by parts formula, we get equality 

E{(A^) xy [Ll 2 ---LlX} = \ a i2*{ L ^ Aai - 2 - 3 U [L° a2 ---L° ar ]°} 

+ i E ^^{{A a ^- 2 ) xy Ll 2 ---Ll i _ i Ll i+i ---Ll r }. 

1=2 

Taking into account identity 

E(XYZ°) = EX E(Y°Z) + EY E(X°Z) + E(X°Y°Z) - E(X°Y°) EZ (4.9) 

and performing elementary transformations (see [6] for details), we get the second 
main inequality 



' 2s-2 



D$ (x, y)<\ (M * {x, y)) ^ + i (u(x, y) * D^) , 
+ \D 2 :t^x,y) + \(D^(x,y),D^) 2s _ 2 

^ (Vtov) * — ) 2S _ 2 + 'P-ff- 1 * (-o) 



where we denoted 



_ (2s + 2)(2 8 + l) 

1^9= — Z 



^ = y 2 v / ^2,(a;,|/) 



and = n^E^i^Os.z)- The initial condition for (4.10) is given by the 

obvious equality D 2 2 \x 1 y) — S xy (4:n 2 )^ 1 . Also we accept that D^J = whenever 
r > 2s and that D ( 2s ] = and D ( 2 J = S sfi . 

Consider a triangular domain of integers 

A = {(s,r) : s > 1, 2 < r < 2s} 

and the numbers U s {x,y), V 2 J(x,y) determined in A by the following system of 
recurrent relations 

U a (x,y) = ±(U*U(x,y)) a _ 1 + ±V™ 1 (x,y), (4.11) 
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and 



where 



V^(x, y)<\(lA* D^(x, y)) ^ + 1 (lA{x, y) * Z>«) ^ 
U'J(x,y) = {2s + 2 f S + 1) U s (x,y), W. = -^M,(t,i), 



(4.12) 



t=i 



D s r ) = n 1 X)"=i ^s^C*!*)) an( i tne initial conditions are given by equalities 
Uo(x,y) = 5 xy and T>f\x,y) = 

The main technical result of the present section is as follows. 

Lemma 4.1. TTie family of numbers {U s (x,y), V s r \x,y), (s,r) € A} exists, is 
uniquely determined by (4-11) and (4-12) and is such that 



sup u£\x,y)<5 xy U^ and sup £>£> (x, y) < X>( 

x,j/=l,...,n a;,y— l,...,n 



r) 



(4.13) 



For any /i > 1/12 i/iere ezisis k > such that for all 1 < s < sq with Sq < nn 2 the 
following inequality holds 



¥>(t) + 



h 



n 2 ( 1 _ T )5/2_ 

moreover, there exists C, 1/24 < C < min{2/i/3, 24} smc/i i/iai inequalities 



(4.14) 



C(3r')! n- 2r ' 
C(3r' + 3)! n 



r(l-r) 



-2r' 



r(l-r) 



-(4r'+5)/2 



if r = 2r' 
if r = 2r' + 1 



(4.15) 



hold for all r and s such that s + 2r + 5 < s . 



The proof of Lemma 4.1 repeats almost literally the proof of Lemmas 2.1, 2.2 
and 2.3 of [6], so we omit the computations and explain only the key points of the 
method. One more reason for this is that the main ingredients of this method will 
be used also in the next section, where we present the principal computations that 
are similar for those of the proof of Lemma 4.1. 

First of all, it follows from relations (4.11) and (4.12) that U2 S (x, y) = for x ^ y. 
Then we can consider the diagonal terms only U2 S (x, x). Moreover, we can introduce 
auxiliary numbers U and V that serve as the upper bounds for sup,,, U{x, x) and 
sup x T>(x, x) and verify relations of the form (4.11) and (4.12) and therefore satisfy 
(4.14) and (4.15). 
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Next, assuming that is of the order o(l), it is easy to deduce from (4.11) 
that the leading contribution to is given by [<£>(t)] s . This generating function is 
determined as a solution of the quadratic equation 

<p(r) = 1 + tV(t)/4 



and therefore is given by equality 

t<p(t) 



i-VT 



(4.16) 



(4.17) 



Assuming that T>i 3 ^ — o(T> [ & z> ), it is not to hard to get from (4.15) with the help of 
(4.17) and the identity (see also relation (7.4) of Section 7) 

1 



S(2) 



(s + l)(2s+ l)m s = 



(1 



\3/2 



(4.18) 



that 



p(2) 



(1+0(1)). 



n 2 (1-t) 2 

Finally, regarding (4.17) and (4.18), it is natural to consider the third derivative 
ip'"(r) as a function that determines the next-in-order corrections for IA S and therefore 
to find the optimal the form of (4.14). 

Now let us return to the relation (4.11) written down for U s . Accepting the 
estimate (4.14) to hold, we see that the first term of the right-hand side of (4.11) can 
be rewritten as 



<P(t) + 



hr 2 



n 2 (l-T) 5 /2 



<p(r) - 1 + 



n 2 (l-r) 5 /2 



+ 



4n 2 



(1-t) 



The negative part of the last expression given by $„ (r) = hr 2 (1 — r) 2 n 2 is exactly 
what we need to compensate the contribution C"(l — t 2 ) -2 n~ 2 to the right-hand side 
of (4.11) that comes from the estimate of T>1_ 1 . It remains to control the last term 
h 2 [r 5 (l — -r)~ 5 ] s (4n 4 ) -1 no to be greater than [$(t)„] s . This is true provided the 
positive ratio k = sups 3 /n 2 is sufficiently small. 

Finally, accepting that f>i 5 ^ = o{V^) 1 we get from (4.12) that 



s-ill 

n 4 



1 



(1-rY 



(l + o(l)). 



This leads to the conclusion that the value 
s(3) 



'[(I' 



v-9/2 



] s can be used to estimate 



T>g . The form of (4.15) for general r is dictated by (4.12) and by the detailed 
analysis of the expressions and constants involved into the computations. 
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5 Correlation function terms 

In the present section we study variables 



p£\*,v)= E 

a+/3+7i + ...+7 r _2=2s 



Qt\*,y)= E l E 

a+^+7i + ...+7 r _2=2s 



and 



T? s \x) = E 

a+^+7i + ...+7 r _2=2s 

that we refer to as to the non-crossing, crossing and diagonal terms, respectively. 

Theorem 5.1. Under conditions of Lemma 4-1, there exists \ < niin{«, 2~ 6 } 
such that inequalities 



e{(^ x (a%l;...l;._ 2 }|, x ^ y 

{(A a )° xy (A< ) y yx L° 7l ...L° yr _,}\, x^y 
E{(A")° X (AP) XX L° 1 ---L° lr _ 2 }\, 



sup P^{x,y)< 

x,y— l,...,n 



C(3r' + 1)! 

C(3r'+4)! 
n 2 ^ 2 



(l-r) 2 -' 



( 1 _ T )2r' + 5/2 



sup Q< r >(a:,j/)< 

s,l/=l,...,n 



C(3r')! 
C(3r'+3)! 

„2r' + l 



(l-r) 2 "-'- 3 / 2 



(1-r) 2 



x— l,...,n 



f C(3r')! 


T 






_(1 


_ r )2r'-3/2 


5 


I C(3r'+3) 


! 


T 




I, n 2 -^ 1 


_(l-r) 2 '-' + 1 


S 



if r = 2r'; 
if r = 2r' + 1, 

if r = 2r'; 
, ifr = 2r' + l, 

if r = 2r'; 
ifr = 2r' + l, 



(5.1) 



(5.2) 



(5.3) 



ftoZd /or aZZ s and r such that s + 2r + 5 < s with Sq < %n 2 . 

We prove Theorem 5.1 by using a modification of the recurrent relations method 
developed in [6]. To derive these relations, let us consider expression 

H(a, (3, r r „ 2 ) = E |(A")° b (A?)° cd ■ ■ ■ L°, <r _ 2 } 

n 

= ^{A at (A^u ((^)° d L;...L; r _ 2 ) } 
t=i 

and apply to the last mathematical the integration by parts formula (7.1). We get 
equality 

H(a,f3,T r - 2 ) = \ ]T>{(.4«- 2 ^) Qb Lj ((A^r cd L; i ---Ll_ 2 )°} 
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+i E E {(^)- (^- 2 - j ) c6 l; • • • l;_ 2 } 
+i X> E {(^ a+7! - 2 u ^ ■ ■ ■ ^ • • • K_ a }. 

Using (4.9), we obtain that H(a 7 (3, r r _ 2 ) can be represented as a sum of ten terms, 

10 

E { (A«)° ab (A?y cd L° 7i ... L° lr _ 2 } = J [k) 0, ft T,- 2 ), (5.4) 

k=i 

where 





3=0 




J(2) : 


= i E (°. & ) E { (^)c<i ^ ^ • • 

3=0 




J(3) 






J(4) = . 


a-2 

--£E{(A«- 2 -^ fe (^)° d } e{l°l; 

j'=o 












= ^E^- 2 -;M) E {( Aj )^ L ° 1 -- 




J^ = 


= ^-E^( a ' d ) E {( Aa+/3 " 2 " j )c6 ^r- 

3=0 




J(8) 


3=0 


z, °.-2 } ; 



J(9) = 4^2 X> ^-Hr,-2(a, 6) E |(A^)° d ^ • • • L^ +1 • • • ^ r _ 2 } , 



J(10) = h X> e {(A«+^- 2 )° b (^)° d l; • • • l;^ . . . ^ r _ 2 } . 

1=1 

Basing on (5.4), we derive a system of recurrent relations for the terms P, Q and 
T. The structure of these relations resembles very much the one of (4.12) and the 
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triangular scheme of recurrent estimates can be applied to them. In general words, 
we assume the estimates of the form (5.1), (5.2) and (5.3) to be true for all the terms 
of the right-hand sides of the relations we get, and prove that the terms of the left- 
hand sides also obey the corresponding estimates. More detailed discussion of the 
triangular scheme of the proof of the recurrent relations can be found in [6] . 

5.1 Relations and estimates for P^l\x,y) 

Taking into account that U2 S {%, y) = when x ^ y and denoting U2 S {x, x) = U2 S (x), 
we get from (5.4) the following system of relations for P^l\x, y); 



2s-2 



+ 



+\p^2-\x,y) + \(pW(x,y)*D^ 
^(u» { x)*D^) 



\(y(x)*DV{yj) 



2s-2 



' j.s-2 An 2 



2s-2 2n 
2 (2s + 2)(2s + l) (r _ 1} 



2s-2 



(x,y) (5.5) 



with the initial condition P^\x, y) = 0. 

Let us consider the case of ■ Regarding the estimate of U2 S {x) (4.6) and 
taking into account formula (7.5) of Section 7, we can write that 



18h 



(1 - t) 3 / 2 



(1 - t) 9 / 2 



(5.6) 



Using these estimates and substituting (5.1) and (5.2) to the right-hand side of (5.5), 
we get the ten terms that present as the following sum; 



where 



p(2r;l) 
2s 



(2r;k) 
2s i 



fc=l 



C(3r + 1)! 
4^ 

lC(3r + l)! 



tV(t) 



(1-r) 

T 



2r 



1 C(3r 

2 n 2 



D! 



(1-r) 



2r-l/2 



v. 



(2r;2) 

2s — 



C(3r)\ 
An 2r 
1 C(3r)! 

9 -n2r 



t 2 <p(t) 



(1-r) 
r 



2r 



(1-r) 



2r 



1 C(3r)! 

2 n 2r 



? (2r;3) = MC(3r + 1)! frC(3r)! C(3r + 4)! 

2s I /l„2r+2 4 n 2r+2 /l„,2r+2 



V 4n 2 ''+2 



4n 2r+2 



(!_ T )2r-l/2_ 

r 4 

_(l_ T )2r+5/2 
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(2r;4) _ ( C 2 4!(3r)! feC(2r-2)(3r-3)i 



2s 



4n 2r+2 



,(2r;5) _ sC( 3r ) ! 



2s 



2n 2 '' 



(1 - r) 21 - 3 / 2 



4n 2r+2 / L(l - r) 2r + 5 / 2 



C(2r-2)(3r-3)! 



+ 



4n 2r 



(1 - r) 2 '- 1 /2 



and 



,(2r; 6 ) _ C(2r - 2)(3r + 1)! (2s + 2)(2s+l) 



2s 



4n 2r+2 



(1 - r) 2l '+ 1 /2 



Regarding the sum of an d ^23" ' 2 \ we can write that 



(2r;l) + ^(2r;2) _ C(3r + 1)! 



2 s 



2 s 



,2r 



where 



An — 



3CV(3r) 



2n 



2r 



(1-r) 



2c 



and Jo = 



C(3r + 2)(3r)! 



■*b - Jo, 



(5.7) 



2n 



2r 



(!_ r )2r-l/2 



The first term of the right-hand side of (5.7) reproduces the expression that we use 

(2r) 

as the estimate of P 2s • Therefore all that we need is to check that the sum of the 
remaining four terms V is less than the sum X +3V In fact, we are going to compare 
the sum of these four terms with Jo- 

Using identity (7.3) (see Section 7), we can write that for any integer k 



,1+k 



(1 - r) 2r+5 /2 



< 



(1 - r) 2r+5 /2 



_ (2s + 4r - 3) (2s + 4r - l)(2s + 4r + 1) 
~ (4r - l)(4r + l)(4r + 3) 

Similar computations show that 



(1 - r) 2 '- 1 / 2 



< 



(1 - T )2r-3/2_ 

Finally, we get equality 



(1 - r) 2 -- 3 / 2 



4r-3 
2s + 4r - 5 



(1 - T ) 2r -V 2 



( 1 _ r )2r+l/2 



2s + 4r - 3 
4r- 1 



(!_ r )2r-l/2 



(5.8) 



(5.9) 



(5.10) 



Using these three relations and taking into account that s + 2r + 1 < s = \r?^ , we 
can write that 



1 ^ r (2r-k) < 16fe X 4 (3r + l)(3r + 3)(3r + 4) 192C X 
Jof^ 2s ~(4r-l) 3 X (4r- l)(4r+ l)(4r + 3) (3r + 2)(4r - l)(4r + 3) 



(4r-3)s 



r - 1 



18h X 



(3r + 2)(2s + 4r-5) (3r - 2)(3r - l)3r(3r + 2) (4r - l)(4r + l)(4r + 3) ' 
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Taking into account that r > 1 , we see that the right-hand side of the last inequality 
is strictly less than 5/6 when h and x verify conditions of Theorem 5.1. The estimate 
(5.1) is proved for P^ ■ 



Let us consider relation (5.5) for P^ s r+1 ^ 
and (5.6) and using expression (4.17), we get an inequality for P 2s 
ten terms that we present as the following sum 



Taking into account estimates (4.6) 
^ that counts 



where 



Pt +1 \*,y)<Y,V: 



(2r+l;fe) 
2s ) 



(5.10) 



fe=i 



■>(2r+l;l) 
2s 



C(3r + 4)! 

4n 2r + 2 
C(3r + 4)! 

2n 2r + 2 



T 2 ip(r) 



(1 - r) 2r+5 /2 



(1 _ r )2r+5/2 



C(3r + 4)! 
2n 2r + 2 



(1-r) 



2r+2 



T)(2r+1;2) 
' 2s 



C(3r + 3)! 

4n 2r + 2 
C(3r + 3)! 

2n 2j -+ 2 



r 2 </>(r) 



(1 - r) 2r+5 /2 

T 

(1 _ T )2r+5/2 



C(3r + 3)! 
2n 2r + 2 



(2r+i;3) _ ( 9Ch{3r + 4) ! ^ 9Cfe(3r + 3)! 9/iC(2r - l)(3r)! 



(! _ T )2r+2_ 
^4 



2s 



^ 2n 2 ''+ 4 + 2n 2r + 4 



2n 2r+4 



(1-r 



|2r+5 



-p(2r+l;4) 
' 2s 



C(3r + 4)! C 2 4!(3r)! 



p (2r+l;5) = sC(3r + 3)! 
2s 2n 2r+2 



4n 2r+2 

r 



4n 2r+2 



(1-r) 



2r+2 



(1 - r) 2l '+! 



C(2r- l)(3r)! 

4n 2r+2 



(1-r) 2 - 



+2 



and 



■P. 



(2r+i;6) _ C(2r - l)(3r + 1)1 (2s + 2)(2 s +l) 
' 2s 2n 2r + 4 ' 2 

Regarding the sum of two first terms, we can write that 



(1-r) 



2r 



where 



and 



- D (2r+l i l) , ^(2r+l;2) _ C(3r + 4)! 
'2s + ' 2s - 2n 2r + 2 

_ C(3r + 5) (3r + 3)! 
2n 2r + 2 

C(3r + 3) (3r + 3)! 
2n 2l -+ 2 



(1 - r) 2r+5 /2 



^i-yi, 



2r+2 



yi 



(1-r) 



(1 _ T )2r+5/2 
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Using identity (7.3) of Section 7, it is not hard to show that for any integer s > 1 



(1 - T ) 2r + 5 / 2 



> 



(1 - T) 2r + 2 



Then we can write that 



x\ + yi < 



C(3r + 3) (3r + 3)! 



,2r+2 



(1 - r) 2r + 2 



Now it remains to show that the sum X)l=3 "P^ 2r+1 ' fe ' does not exceed A^. To do 
this, we will use (7.6) and its consequences given by the following three relations; 



(1-r) 



2r+5 



< 



(2r + 2)(2r + 3)(2r + 4) 
2r + l " 



(1-r) 



2r+2 



(1 - Tp+l 



and 



(1-r) 



2r 



s + 2r 
2r(2r + 1) 



(1-r) 



2r+2 



(s + 2r- l)(s + 2r) 



(1-r) 



2r+2 



(5.11) 
(5.12) 

(5.13) 



Using these three relations, we can write that 



1 



(2r+l;fe) 



2s 



< 



9hx(3r + 5) 



fe=3 



3C 



4(3r + 3)(2r + 2)(2r + 3)(2r + 4) 



2r + 1 

+ 7777, 777 + 



4(3r + 5) (3r + l)(3r + 2)(3r + 3) 2 2(3r + 5) (3r + l)(3r + 2)(3r + 3) 2 

9/ix(2r-l) r(2r- l)(2r+ 1) 

+ 2(3r + l)(3r + 2)(3r + 3) 2 + (3r + 2)(3r + 3) 2 ' 

It is clear that under conditions of Theorem 5.1 the right-hand side of the last in- 
equality is strictly less than 1 for any r > 1. The estimate (5.1) is proved for the 
variable P 2 ^ r+1 '(x, y). 

5.2 Relations and estimates for Q^ix^y) 

Regarding (5.4) with a = d = x and b = c = y and taking into account equality 
U2s{x,y) = for x 7^ y, we obtain the following recurrence; 

Qi\x,y) < \ (M * Q«(x, y )) 2s _ 2 + Jq^W) 



2s-2 



2s-2 



16 



s - w ^ ■ (r ~ 2) (25+2)(25+1) Q^r i 2 ) (x, y ), 



(5.14) 



where we denoted 

[4 (a;) = (2s + 2)(7 2s . 

Let us consider the case of ■ Regarding (4.6) and using formulas (7.4), we 
obtain that 

1 hr 2 

+ 



u'M < 2 



vT^T n 2 (l-r)7/5 



(5.15) 



Using this estimate and substituting (5.1) and (5.2) to the right-hand side of (5.14), 
we get twelve terms that we regroup in the sum of six terms as follows; 



Qr/W)<£a 



(2r;fc) 
2s i 



fc=l 



where 



(2,;i) C(3r)! 

C(3r)! 



rV(r) 



2n 2r - 1 



(1 - Tp- 3 / 2 

T 



(1 - T )*r-3/2 



C(3r)\ 



(1-r) 



2r-2 



(5.16) 



^2s 



c 



4n 2r+l 



(/i(3r)! + (3r + 3)! + C3!(3r)! + 2h(3r - 3)!) 



D (2r;3) = C(3r)! 
^2s 2n 2r+l 

( 2r: 4) _ Ch 
^ 2s ~ 2n 2r + 3 



Q 



(l-rr 

*(3r-3)! (1 _ r)2r+4 
(2r;5) _ C(3r - 3)! 



+ (r-l)(s + l)(2s-l) 
-s(3r)! 



_(l_ T )2r+l 

r 



(1 - r) 2 '- 1 

^4 



(1-r) 



2r+3 



2s 



2n 2r - 



and 



Qr 6) = ^T(^,! 



Regarding Q^' X \ we can write that 



(1 - r) 2l '+ 1 /2 



(1 - t) 2 '-3/ 2 

+ /i(3r-3)! 



(1 - r) 2r + 3 / 2 



n (2r;l) _ g(3r)! 



,2r-l 



(1 - r) 2 '- 3 / 2 



* 2 - y 2 , 



where 



C(3r)! 



.(1-r)* 



^2 



g(3r) 
2n 2 '- 



_(l_ T )2r-3/2_ 



(5.17) 
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We are going to show that the sum J2k=2 3 4 Q^s'^ IS strictly less than X 2 and that 

the sum 0,2s ~\~ 0,2s is strictly less than 3^2- To do this, we will use the following 
consequences of formulas (7.3) and (7.6): 

_ (s + 2r- 3)(s + 2r - 2) • • • (s + 2r - 4 + q) 



(1 - T ) 2r - 2 +i_ 
for q = 1,2, 3, 5, 6 and 



(1 - r) 2l '+ 1 /2 



(2r-l)2r---(2r-2 + 9 ) 
(2s + 4r-5)(2s + 4r-3) 



(1-r) 



2r-2 



(4r-3)(4r-l) 



(1 - T ) 2f -3/ 2 



(5.18) 



and 



(1 - r) 2r + 3 / 2 



(2s + 4r - 5)(2s + 4r - 3) (2s + 4r - 1) 



(4r-3)(4r-l)(4r + l) 



(1 - r) 2 '- 3 / 2 



Using (5.17); we conclude that 



(5.19) 



(2r;fc) ^ (3r + l)(3r + 2)(3r + 3) + h + 3C 



fe=2,3,4 



+ - 



2r(2r- l)(2r + 1) 

fax 2 fr 2 x 2 

7! 8! ' 



2r(2r - 1) ' A ' 7! ' 8! ' (5 - 2 ° } 

where the right-hand side is strictly less than 1 for any s, r > 1 and sj^/n 2 < Also 
we have inequality 



UQf> h) + Qf fi) )< 



1 



y 2 



+- 



3r(3r- l)(3r-2) (4r - l)(4r - 3) 
hx 



(4r - l)(4r + l)(3r)(3r - l)(3r - 2) ' 
where the right-hand side is strictly less than 1 under conditions of Theorem 5.1. 

Let consider Q^" +1 \x, y). Wc have from (5.14) that 



(5.21) 



where 



(2r+l;l) 
«2s 



fe=i 



C(3r + 3)! 

4n 2r+l 

C(3r + 3)! 
2n 2r +! 



tV(t) 



(1 - r) 2 ''+! 

T 

(1 - T) 2 ^ 1 



C(3r + 3)! 
2n 2r + 1 



(5.22) 



(l_ T )2r+l/2 
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C 



-((3r + 3)! + 6C(3r)!/n) 



+C(3r)! 



4n 2r+l 

(2r-l)(« + l)(2« + l) 



(1 - r) 2r +V2 



4n 2r+l 



2^ +1;3) = ^(M3,)! + (3r + 3)!) 



(1 - t) 2 ^ 3 / 2 

T 



- ^3(C(3, + 3)! + (3r)!) 



(1 - T ) 2r +V 2 



(1 - T ) 2r + 5 / 2 
h 2 C(3r)\ 



+ 



2n 2r+A 



(1 - r ) 2l '+ 13 / 2 



and 



(2r+l;5) = C(3r)! 
^2s 2n 2r+l 



(1 - r) 2 ^ 1 



(1-r) 



2r+4 



Regarding Q^ +1,1 \ we can write that 



n (2r+l;l)_ C(3r + 3)! 
^2s 



,2r+l 



where 



X 3 



C(3r + 3)! 
2n 2l '+! 



(1 - r) 2 ^ 1 



(1 - r) 2 ^ 1 

C(3r + 3)! 



x 3 - y 3 , 



(5.23) 



y 3 



2n 2r+1 



(1 - r) 2r + 1 /2 



Now it is sufficient to show that J2k=2 3 4 Q^s^'^ i s strictly less than 3^3 and that 

(2r+l;5) 



2.s 



is strictly less than X 3 . 



It follows from (7.6) that 

T 



(1 - r) 2 '- 3 / 2 



(4r-3)(4r- 1) 



and that 



(2s + 4r- 3)(2s + 4r - 1) 
(2s ) 9 



(1 - r) 2l '+ 1 /2 



(X _ T )2r+l/2+q 

Then we easily get inequality 



< 



(4 r )s 



y 3 



(2r+l;fc) ^18 

2s S 2 27 



_(! _ r )2r+l/2_ 



16 



49 ' 



(5.24) 



fc=2,3,4 

Clearly, the right-hand side of it is less than 1 under conditions of Theorem 5.1. 



Using (7.6), we can write that 
1 Q (^)< 



2rs 



2 h\ 

+ T7-. — To + 



X 3 " (s + 2r)(3r) 3 ' (3r) 3 ' (2r) 3 ' 



(5.25) 



The right-hand side of this inequality is obviously less than 1 under conditions of 
Theorem 5.1. Inequality (5.2) is proved. 
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5.3 Relations and estimates for t![ s \x) 

Using (5.4) with a = b = c = d = x, we get the following recurrence for the diagonal 
termT 2 ( s r) (z); 



1 2s 



(x) <\(M* tW(x)) 2s _ 2 + \ (U(x) * D {r) (x)) 2s _ 2 + \t£?(x) 
T£U X )*D^) 2s _ 2 + l(u>(x)*U { x)*D^) 2s _ 2 



(u { x)*D^ { x)) 2s _ 2 + ^T£_ 2 (x) 
D (r-2)\ + r_l . (2s + 2)(2s + l) r(r _ 1} 



The right-hand side of (5.26) contains more terms than those of the relations for the 

(r) (r) 

non-crossing terms P 2s ' and the crossing term Q 2s . Therefore the total number of 
terms to consider raises up to 16. However, the estimates repeat in the most part 
the estimates performed to prove (5.2). This is because the second and the eights 
terms of the right-hand side of (5.21) that are absent in (5.14) but present in (5.5) 
are of the order smaller than the leading terms of the right-hand side of (5.21). That 
is why the diagonal term T 2 J{x) is bounded by the same expression (5.3) as the 
non-crossing term Q 2 ] (5.2). 

We do not present the detailed proof of (5.3) because it is very similar to that of 
the proof of (5.2) and uses the same formulas of Section 7. Indeed, when estimating 
T 2 ^ r \ we conclude from relation (5.25) and expressions (5.3) that the leading term 
and the negative part —X\ — 3^4 are given by the corresponding terms of the right- 
hand side of relation (5.17). Then it is not hard to see that the "extra" terms coming 
from the right-hand side of (5.26) with respect to (5.14) add the terms 

X . h 2 X 2 , hx , h X 3 

— + — — — — — and —— + 



24 3r(3r-l)(3r-2) 3 12(3r - l)(3r - 2) 

to the right-hand sides of (5.20) and (5.21), respectively. Certainly, this does not 
alter much the result of the sum that is still strictly less than 1 under conditions of 
Theorem 5.1. 

Regarding T 2s , we see that the leading term and the negative contributions 
A5 and 3^5 are exactly the same as the corresponding terms of the right-hand side of 
(5.23). The "extra" terms of (5.26) produce then the additional terms 

% and £ + * + 



(4r+l)(4r + 3) 32 8(3r) 3 8(3r) 3 

for the right-hand sides of (5.24) and (5.25), respectively. Certainly, this does not 
alter the result of sums that are strictly less than 1 under conditions of Theorem 5.1. 
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6 Proof of Theorem 2.1 



Regarding variable R n {s',s") (4.3), we apply the result (4.6) of Theorem 4.1 to the 
factors of the first term Rn\s', s") and write that 



R£\-s',s")<n> 



^(r)+hT 2 n- 2 (l-T)- 5 ^ ■ ^(r)+hr 2 n- 2 (l-T)- 5 ^ 

J s' L 

where s = s — 1. Repeating computations of (4.19), we get inequality 
R£\s',s") < 16n 2 L(T) + hT 2 n- 2 (l-r)- 5 / 2 ]_ ■ Ur) + hr 2 n - 2 {\ - r)" 5 / 2 



his') 

< 16n JTigi m s n [ 1 + 



/ his 1 ') 3 
1 + 



It follows from the expression for m s (4.5) that m s — (7rs 3 ) _1 / 2 (l +o(l)) as s — > oo. 
Then in the limit (2.5) we have the bound 



li n,sup fl y ( ,-,,-0<16 " + ^^r 

rwoo TTVX X 



(6.1) 



Let us pass to the last term of the sum (4.3). Using the results of Theorem 5.1, 
we can write that 

n n 
x,y=l 



x=l 



24s'C\ /24s"C 



6s'C 



+ n m s '_i 



6s'C 



m s »_i 



J \ J \ n / \ n 

Taking into account the asymptotic expression fro m s , we get in the limit (2.5) 

*£V,*")< {^) (xV0 1/3 (l + O (l)) + n5/3 ^ /)1/6 (l + o(l)). (6.2) 



Similar computations show that 



WO + «.?»(»'.»») <- % (»±W + <*2$£M) . ( o„ 

Remembering the factor V^/n 2 of (4.3), we see that S„(s',s") (3.3) vanishes in 
the limit (2.5) as n — > oo. 

Let us consider the variable S n (s',s") (3.7) and its representation in four terms 
similar to (4.3). It follows from the results of Section 4 that the terms S^\s' , s"), 
fc = 1,2,3 that contain factors U a (x,y)Up(x,y) are equal to zero. Then only the 
term 



W,*") = ]T £ E{(^)°(A*)°} ]T E{(^)°(A*)°}), 

x,j/=l \qi+^i=2s' / \a 2 +/3 2 =2s" 
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gives a non-zero contribution to S n (s / ,s / '). Repeating computations of Subsection 
5.2, it is not hard to show that «sl 4 ^(s', s") = o(l) in the limit (2.5) as n — >• oo. This 
completes the proof of Theorem 2.1. 



7 Auxiliary relations 

For completeness, let us refer to some of the equalities and identities of [6] that we 
use in the present paper. 



The first equality is a consequence of the integration by parts formula applied 
to the normal (Gaussian) random variable £ ~ j\f(0,v 2 ) that is E£/(£) = w 2 E/'(£), 
where f(x) is a non-random function such that corresponding mathematical expecta- 
tions exist. Then for the random matrix A from GUE we get, in particular, relation 



fe-i 



(A k ) uv = ^ X>{(^U (A k -^) xv ) . 



(7.1) 



3=0 



The second identity relates the generating functions (1 — r) r x l 2 with the mo- 
ments m s . It is not hard to obtain that 



1 



(1 - ry+v* 



(2r+2s\ 



(ZD 



or in equivalent form, 



(1 - t)''+ 1 /2 



1 (2r + 2s)! r! 



2 2s s! (r + s)\ (2r)\ 



Two particular cases are important: 

(2s + l)(2s + 2) 



■ m s = 



and 



(2s + l)(2s + 2)(2s + 3) 



(1 - t) 3 / 2 
1 



We also use the equality 



1 



(1-r) 



3! 



fc+i 



■ TO, 



_(l-r) 5 /2 
(s + l)(.s + 2)---(s + fc) 



(7.2) 



(7.3) 



(7.4) 



(7.5) 



(7.6) 
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